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Abstract 

We demonstrate using the classical density functional theory (DFT) model that an intermediate 
degree of mass localization in the amorphous state is essential for producing the boson peak. The 
localization length I is identified from the width of the gaussian density profile in terms of which 
the inhomogeneous density n(x) of the solid is expressed in DFT. At a fixed average density, there 
exists a limiting value £q of £ signifying a minimum mass localization in the amorphous state. For 
more delocalized states {£ > £q) occurrence of boson peak is unfeasible . 







Studies of inelastic scattering of light and neutron from amorphous solids show that 
the disordered systems are characterized by an excess density of vibrational modes over that 
predicted by the standard Debye distribution. This excess in the density of states (DOS) g(uj) 
for a disordered system typically occurs at the THz frequency (u) range and is identified as 
the so called boson peak in the q(u)/u) 2 vs. u> plot. The boson peak height decreases and its 
location shifts towards the higher frequencies with the increase of pressure or density PJ-[2[. 
A number of different theoretical models have been developed to explain the anomalous low- 
energy excitations [^-[h]]. It is now generally believed that the boson peak is a manifestation 
of disorder and is key to the understanding of the vibrational states of glassy materials. The 
origin of boson peak in the disordered solid has been often linked to transverse sound modes 



for the system 



111 ]. It has been argued that the boson peak results from a band of random 



transverse acoustic vibrational states [7|. Extensive computer simulations of several glass 
forming systems 12| also indicated subsequently that the origin of the boson peak is the 
transverse vibrational modes associated with defective soft structures in the disordered state. 
In a recent work, Chumakov et al. |l3] explain the boson peak in terms of (predominantly) 
transverse sound waves. 

A common aspect of a class of the boson peak models mentioned above is that they 
are based on the existence of localized modes in the amorphous solid [14]. Longitudinal 
acoustic phonon modes are observed in a liquid by ultrasonic, optical, or inelastic scattering 
experiments. These hydrodynamic modes occur due to the conservation laws for the system. 
For the crystalline state the isotropic symmetry of the liquid is broken and the constituent 
particles vibrate around the sites of a lattice with l ong range order. This symmetry breaking 
leads to the development of the Goldstone modes [l5l] or the transverse sound modes in the 
crystal in addition to the longitudinal sound modes. For the amorphous solid the symmetry 
is not broken over long length scales like that in a crystal. When the atomic vibration 
wavelength in a liquid approaches the atomic nearest neighbor distance, i.e., in the tera 
hertz frequency region, there is a solid like cage formation giving rise to a restoring force for 
the acoustic transverse modes. 

A typical choice of the order parameter for the crystalline as well as the amorphous state 
is the inhomogeneous density n(r) with smooth spatial dependence. The density for the 
solid state is expressed in terms of gaussian profiles of width i and centered on the sites of 



a chosen lattice {R,:} 



17|. 
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with 1 / yja = £. The coarse grained picture for the solid presented above involves two inher- 
ent length scales. First, the microscopic length scale a (say) associated with the interaction 
potential between the constituent particles. For a hard sphere potential, o corresponds to 
the hard sphere diameter. Second, the length £ which represents the degree of localization 
of mass in the system. The homogeneous liquid state is characterized by the limit aa 2 — > 
or £ » a. The sharply localized density profiles of the crystalline state correspond to the 
limit aa 2 » 1 or £ << a. The ratio £ja is generally referred to as the Lindeman parameter 
for the crystalline. In the present work we compare the free energy contribution due to the 
vibrational modes in the amorphous state to the corresponding density functional result for 
the ideal gas part /id[n][16] of the free energy for an inhomogeneous density n(r). We show 
that for a a chosen average density, the characteristic gaussian profile of n(r) signifies a 
minimum localization length £q below which the appearance of boson peak in the reduced 
density of states is unfeasible. 

The equilibrium state of the fluid is identified by minimizing a proper thermodynamic 
potential or free energy with respect to the density n(x). The free energy is expressed as 
a sum of two parts which are respectively the free energy of the non-interacting system 
and the contribution F ex due to interactions between the particles. Fa is obtained from the 
logarithm of the partition function (V / A 3 ) N / N\ for the noninteracting system of N particle, 
where A = fij ^iTrmk^F denotes the thermal wave length at temperature T. Thus $F\& = 
Vn Q (ln(n A 3 ) — 1) for the uniform density fluid n . For the inhomogeneous density n(x), 
the expression for the free energy is easily generalized to = / dm(r)[ln{n(r)A 3 } — 1]. 

For the uniform density n changing to a nonuniform one n(x), the ideal gas contribution 
for the N particles system changes by an amount n j drn(r)[lnn(r)], where n(r) = n(r)/no- 
Since by definition n(r) is always positive, using the Gibbs inequality {xlnx — x + 1} > for 
positive x, it follows that entropy drops upon localization of the particles. This is a result 
of the restriction of available phase space. If a is larger than a limiting value «o ( sa y) the 
density n(r) is well approximated in terms of the gaussian profile centered at the nearest 
lattice site. The free energy per particle fid for large aa 2 is obtained as, 
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/idN«4 + 31n (Vf A )- (2) 

In a microscopic approach, vibrational modes with localized density profiles contribute to 
the free energy and the latter is obtained in terms of the density of vibrational states g{oj) of 
the system. The free energy of the non-interacting system is obtained as F ld = — /ceTlnZ/v 
where is the partition function for 3N harmonic oscillators in constant NVT ensemble. 
Going from discrete to continuum in the frequency spectrum, we obtain the ideal gas free 
energy per particle in units of (3~ l as /3F id /N = K,(x m x)g(x)dx = f [g, x m (a)], where 
k(x) = —\{x + ln[l — e~ x ] — 2x/(e x — 1)}. u m is the upper cutoff of frequency up to 
which g(oo) is nonzero and it is expressed in terms of the dimensionless quantity j3hu m = 
x m =x m (a). The scaled density of states g{x) is obtained from g(uo) as g{x) = (3n ) 1 u m g(uj) 
where x = u/u m is the reduced frequency. The expression for k(x) indicates the functional 
dependence of the integral fo on the scaled density of states g. There are 3uq vibrational 
modes for an average no number of particles in a unit volume giving rise to the normalization 
condition J^ m g(u)du = 3n for g(u). The cutoff frequency u m is a characteristic property 
of the material and sets the shortest length scale for the vibrational modes. For the Debye 
distribution, g{u) = <7_d(w) = (9n /cjf))w 2 is nonzero in the range < uj < Ue>. The 

1 /3 

Debye frequency up is obtained as wo = (67m ) ' c, with c being the speed of sound. The 
scaled (Debye) distribution g D as gn(x) = 3x 2 where x = cu/up. Considering the case 
x m =XD = (31tjjJd « 1? analysis of the integral fo[g,x m (a)] for g(x) = go(x) obtains the 
following result to leading orders in x^: 
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fo[g D ,x D ] = -- + 3\nx D + 0(x 2 D ). (3) 

This mapping identifies the Debye cutoff frequency uo in terms of the corresponding width 
parameter a of density functional description as (3hujD=XD(a) = \/ a/nA. For a given A, a 
corresponding a is identified such that for all a > «o the asymptotic form (j2j) will represent 
the ideal gas free energy. The value of A/cr should be such that xd(«o) < where the 
asymptotic form (j3J) for fo(x) holds for x < xl- 

In DFT the free energy /id (a) for the amorphous state is evaluated numerically using 
the n(r) corresponding to the parameter value a. The centers for the gaussian density 
profiles {Rj} are distributed on a random lattice and the corresponding site-site correlation 
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function g(R) is obtained from the Bernal's [18| random structure generated through the 
Bennett algorithm jl9fl . We approximate the g(R) through the following relation g(R) = 
gB[R(Tj/rj ) lj/3 ] where r\ denotes the average packing fraction and r] is used as a scaling 



parameter for the structure 



20] such that at rj = rjo Bernal's structure is obtained. Smaller 



the parameter a in n(x), larger number of sites are needed to be included for an accurate 
evaluation of the f 1( i- The free energy evaluated numerically agrees with the corresponding 
asymptotic expression (jSJ) for a > ao- This is shown in the inset of Fig. 1. The difference 
A/id between the two evaluations of /m is shown in Fig. 1 for a < ao range. This result 
corresponds to hard sphere system with the fixed average density n a 3 = 1.1 or packing 
fraction (p = 7m o" 3 /6 = .576. The thermal wavelength is kept constant at A/a = 0.0137. 
The free energy obtained from the DFT expression /^[^(r)] corresponding to n(r) in the 
range a < ao, deviates completely from the asymptotic result (J2J). It extends to the correct 
a— >0 limit (ln{n A 3 } — 1) for the uniform liquid state. A key observation here is that 
the free energy curve for small a < ao matches with the corresponding result obtained 
from the microscopic expression fo[g,x m (a)] with a modified density of states which is 
different from the Debye distribution. The normalization condition for g(uj) is maintained 
with a corresponding uj m {a) or equivalently x m (a) different from that for the Debye case (= 
\Ja~pn/\). Therefore we require that the function x m (a) — » xd(ck) as a — > ao- The correction 
part in this functional relation is denoted by C(a), such that x m (a) = i£)(a ) + C(a) with 
C(a)— >0 as a— >ao- Similarly, for the scaled distribution function g(x) we use a correction 
over the corresponding Debye form go{x) = 3x 2 in terms of the function A, such that 
g(x, x m (a)) = 3x 2 (l + A(x,a)). Since for all frequencies x, we must have g — > g~z> as 
a— >ao, we express A with the separation of variables. Thus A(x,a) = B(a)A(x) where 
B(a)— 7-0 as a— >ao- Furthermore in order to maintain the normalization condition of g(uj), 
the function A(x) for all a must satisfy x 2 A(x)dx = 0. The two functions C(a) and 
B(a) are parameterized in terms of two respective polynomials of the separation parameter 
e = 1 — a/ao. These functions are determined by requiring that the fo[g,x m (a)] obtained 
with the corresponding g and x m agree with that of the DFT expression for in the range 
a<ao- We choose the x dependence of the A(x) so as to have an intermediate peak over the 
whole frequency range. The relative frequency x at the peak of A(x) is a parameter which 
is kept fixed through out this work. In order to maintain the positivity of the density of 
states the function A (a;) is also kept confined in the range ±1. In Fig. 2 we show the results 
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obtained for the appropriate density of states g(x) which satisfy the above constraints of 
density functional theory. The reduced density of states g(u)/uj 2 vs. u/ojd is shown in Fig. 
2 for three different values of the localization parameter t = 1/y/a. The height Kb of the 
boson peak increases as i grows and is shown in the inset of Fig. 2. Larger i or smaller a 
corresponds to increased spreading of the density profiles and hence decreased localization. 

As the density profiles get increasingly spread out, i.e., as a decreases, the intensity or 
height of the boson peak continues to grow while upper cutoff u m of the frequency keeps 
decreasing. However, this trend does not continue indefinitely towards a— >0. To understand 
this behavior we note that the upper cutoff u m sets the shortest wavelength A m i n for the 
sound waves. In a valid continuum description of the solid, A m i n cannot get smaller than 
the average mean square displacement of the individual particles. In the density functional 
theory formulated in terms of gaussian density profiles, the average particle displacement 
is £ = l/\fa and this is a reasonable estimate for A m i n . The corresponding upper limit in 
the frequency of the sound wave is obtained from the relation fl L / (2nc) = yfa. where c is 
the speed of sound. In Fig. 3, for n§o z = 1.1 this limiting value for the upper limit vs. 
cutoff \J~OL is shown as a dashed line. The corresponding upper cutoff of the frequency w m (a) 
obtained by matching microscopic expression for the free energy with the density functional 
expression (as described with respect to Fig. 2) is shown as a solid line. We mark with 
an arrow in Fig. 3, the point a = a m i n where the two lines cross. This point depicts the 
maximum delocalized state for which the boson peak in g(u) is observed in the vibrational 
density of states. For more delocalized states the description in terms of sound waves breaks 
down. In the inset of Fig. 3, we show how the point of maximum derealization obtained in 
the main figure changes with the density of a hard sphere system. 



Metastable states of strongly heterogeneous density distribution has been obtained 
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23| with free energy minimum intermediate between the uniform liquid and crystals with 



long range order. We regulate the parameter rj Q to define specific sets of random structures. 
The low a or partially delocalized state exists only if total free energy is computed with the 
n(x) for lattice points {Ri} on a random structure. Thus disorder is essential for locating 
the metastable minimum of the free energy in the delocalized region (a < ao). This is due to 
the interaction contribution to the free energy. However the corresponding F^ in the small 
a region remains the same with {Rj} representing a regular crystalline lattice with long 
range order instead of the Bernel structure. With = exp[— Kf/(4a)} the definition (TjQ) 
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for the density n(r) reduces to the form n(r) = tiq + mje lK ^ r in terms of an expansionf 
involving reciprocal lattice vector (RLV) {K,;} . It is useful to note at this point that the 
boson peak in disordered system has been studied in terms of a geometrically perfect crystal 
having random interactions between the neighbors {?], Q] • 

For the hard sphere system with n^a 3 = 1.1 we obtain the interaction contribution to the 
free energy in terms of standard Ramakrishnan-Yussouff functional [2^| using the solution 



of the Percus-Yevick equation with Verlet-Weiss[27|]. For the heterogeneous glassy state, the 
free energy barrier fg is obtained as the height of the maximum in the corresponding free 
energy curve from the metastable minimum at small nonzero a (inset (a) of Fig. 4). For the 
optimum value of a signifying a metastable state, g(uj) is obtained. In Fig. 4 we show the 
variation of fg for the parameter range .62 < rj < .68, with the corresponding boson peak 
height hg. As the free energy barrier fg becoming weaker, the height of the boson peak hg 
goes down. In the inset (b) of Fig. 4 is shown the fg vs. i = 1/ yfa curve. As the mass 
localization gets weaker, i.e. with decreasing (increasing) of a(£) the barrier height fg gets 
smaller. A lower barrier height implies that the structural degradation is less hindered in the 
corresponding system, i.e., represents a glassy material of higher fragility. The dependence 
of the peak height hg and the peak frequency u p on the pressure are respectively shown in 
Fig. 5 and its inset and conforms to experimental findings. 

An important characteristic of the boson peak is that it gets weaker, more fragile the 
system is. This link between the long-time relaxation behavior (fragility) of the glass forming 
material and the vibrational modes in the THz region, follows in a natural manner in the 
present model. Since less fragile or stronger liquids tend to form network structures, the 
density profiles are more localized for them. Therefore increase of fragility is synonymous 
with decrease of a, i.e., less localized density profiles and lower barriers. Since as discussed 
above the boson peak height decreases with decreasing a, it also implies weaker boson peak 



for more fragile systems 28 



291 ] . Our present approach of DFT is similar to the models of 



disordered solids in terms of springs. For large values of a, the sharply localized density 
profiles are interpreted in terms of harmonic oscillators with spring constant k which is 



proportional to the width parameter a 



30j | . A natural extension of the present model will 



be to include fluctuations of the widths a at different sites on the amorphous structure 
as an appropriate description of the heterogeneous glassy state. PK acknowledges CSIR, 
India for financial support. SPD acknowledges BRNS, India for financial support under 
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FIG. 1: Inset : ideal gas part of free energy /id per particle from numerical evaluation of the density 
functional expression (solid line) and from the asymptotic formula ([2]) (dashed line) . The difference 
between the the asymptotic formula and the numerical result for fa is shown in an enlarged form 
in the main figu 
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FIG. 2: The reduced density of states g(u)/uj 2 in units of c 3 ( c is the sound speed) vs. uj/ujd 
corresponding to aa 2 : 13 (solid), 17(dashed), 21 (dot-dashed), for a hard sphere system at density 
of poa 3 = 1.1. Inset shows corresponding boson peak height Kb vs. localization length I (= 1/yfa). 
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FIG. 3: For a hard sphere system at density poa 3 = 1.1, the ratio uj m juj£, (see text) vs. ^fao (filled 
dots). Dashed line shows the upper limit of the frequency of the sound waves for localized the 
vibrational modes (see text). a m ; n shown with an arrow is the minimum possible value of oj m /oJu 
for the boson peak. Inset shows variation of £o (= 1/ y/a m i n in units of a ) with density po<J 3 for 
the hard sphere system. 
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FIG. 4: Boson peak height Kb vs. the free energy barrier fg for the heterogeneous metastable 
states obtained in the ao > a region corresponding to different structures of a hard sphere system 
of poa 3 = 1.1. Inset a) definition of /b and b) /b vs. the width t of the gaussian density profiles. 
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FIG. 5: Ratio of frequency u p at the boson peak to the Debye frequency ud vs. pressure P of the 
system ( in units of (/3cr 3 ) Inset shows height Tig of the peak vs. pressure P. 
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